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Abstract. We describe and rigorously justify the nonlinear interaction of 
highly oscillatory waves in nonlinear Schrodinger equations, posed on Eu- 
clidean space or on the torus. Our scaling corresponds to a weakly nonlinear 
regime where the nonlinearity affects the leading order amplitude of the solu- 
tion, but does not alter the rapid oscillations. We consider initial states which 
are superpositions of slowly modulated plane waves, and use the framework 
of Wiener algebras. A detailed analysis of the corresponding nonlinear wave 
mixing phenomena is given, including a geometric interpretation on the res- 
onance structure for cubic nonlinearities. As an application, we recover and 
extend some instability results for the nonlinear Schrodinger equation on the 
torus in negative order Sobolev spaces. 



1. Introduction 

1.1. Physical motivation. The (cubic) nonlinear Schrodinger equation (NLS) 
(1.1) idtu + ^Au = \\u\'^u, 

with A G R*, is one of the most important models in nonhnear science. It describes 
a large number of physical phenomena in nonlinear optics, quantum superfluids, 
plasma physics or water waves, see e.g. [30] for a general overview. Independent 
of its physical context one should think of (1.1) as a description of nonlinear waves 
propagating in a dispersive medium. In the present work we are interested in 
describing the possible resonant interactions of such waves, often referred to as wave 
mixing. The study of this nonlinear phenomena is of significant mathematical and 
physical interest: for example, in the context of fiber optics, where (1.1) describes 
the time-evolution of the (complex-valued) electric field amplitude of an optical 
pulse, it is known that the dominant nonlinear process limiting the information 
capacity of each individual channel is given by four- wave mixing, cf. [Ki, 32]. Due 
to its cubic nonlinearity. (1.1) seems to be a natural candidate for the investigation 
of this particular wave mixing phenomena. Similarly, four wave mixing appears 
in the context of plasma physics where NLS type models are used to describe 
the propagation of Alfven waves [2.s]. Moreover, recent physical experiments have 
shown the possibility of matter- wave mixing in Bose-Einstein condensates [12]. A 
formal theoretical treatment, based on the Gross-Pitaevskii equation {i.e. a cubic 
NLS describing the condensate wave function in a mean- field limit), can be found 
in [31, 17]. Finally, we also want to mention the closely related studies on so-called 
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auto-resonant solutions of NLS given in [I'A] where again wave mixing phenomena 
are used as a method of excitation and control of multi-phase waves. 

Due to the high complexity of the problem most of the aforementioned works are 
restricted to the study of small amplitude waves, representing, in some sense, the 
lowest order nonlinear effects in systems which can approximately be described by a 
linear superposition of waves. In addition a slowly varying amplitude approximation 
is usually deployed. By doing so one restricts himself to resonance phenomena 
which are adiabatically stable over large space- and time-scales. We shall follow 
this approach by introducing a small parameter < £ 1, which represents the 
microscopic/macroscopic scale ratio, and consider a rescaled version of (1.1): 

=■2 

(1.2) iedtu" + —Au' ^ Xe\u^^u^. 

This is a semi- classically scaled NLS [G] representing the time evolution of the 
wave field u'^(t,x) on macroscopic length- and time-scales. In the following we 
seek an asymptotic description of as e ^ on space/time-intervals, which are 
independent of e. Note that due to the small parameter e in front of the nonlinearity, 
we consider a weakly nonlinear regime. This means that the nonlinearity does not 
affect the geometry of the propagation, see §1.2 below. Technically, it does not 
show up in the eikonal equation, but only in the transport equations determining 
the modulation of the leading order amplitudes. In view of these remarks, the sign 
of A (focusing or defocusing nonlinearity) turns out to be irrelevant. 

1.2. A general formal computation. In order to describe the appearance of 
the wave mixing in solutions to (1.6), we follow the Wentzel-Krammers-Brillouin 
(WKB) approach, as first rigorously settled by Lax [2.3]. Consider approximate 
solutions of (1.2) in the form of high-frequency wave packets, such as 

(1.3) a(t,x)e"^(*'")/^ 

For such a single mode to be an approximate solution, it is necessary that the rapid 
oscillations are carried by a phase (j) which solves the eikonal equation (see [G] , where 
also other regimes, in terms of the size of the coupling constant, are discussed): 

(1.4) at0+i|V0|2 =0. 

Nonlinear interactions of high frequency waves are then found by considering su- 
perpositions of wave packets (1.3). By the cubic interaction, three phases (pi, (j)2 
and 03 generate 

<^ = - <^2 + <^3- 

The corresponding term is relevant at leading order if and only if this new phase (j) 
is characteristic, i.e. solves the eikonal equation (1.4) while also each j = 1, 2, 3 
does so. More generally, we will have to construct a set of phases for some 

index set J G such that each (f>j is characteristic, and the set is stable under the 
nonlinear interaction. That is, if k,i,m G J are such that (f> = (j)k — (jje ~^ 4>m is 
characteristic, then e {<t'j}ieJ- Given some index j g J, the set of (four-wave) 
resonances leading to the phase is then 

= {{k, £, m) e ; (j)k - cj)t + = 't>j]- 

One of the tasks of this work to study the structure of Ij . A first important step is 
obtained by plugging = 0^ — (/)£ + (/),„ into (1.4), since then, an easy calculation 
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shows that 4> is characteristic if and only if the foUowing resonance condition is 
satisfied: 

(1.5) (V0£ - V(/)™) • (V0£ - V0fc) = 0. 

Obviously this is a quite severe restriction in one spatial dimension, while in higher 
dimensions there are many possibilities to satisfy (1.5). In order to gain more 
insight we shall restrict ourselves from now on to the case of plane waves {i.e. 
linear phases, see §2.1), This choice allows for a more detailed mathematical study 
and is also the most important case from the physical point of view, cf. [31, 13]. 
The precise mathematical setting is then as follows. 

1.3. Basic mathematical setting and outline. In the following the space vari- 
able X E M will either belong to the whole Euclidean space M = M'', or to the torus 
M=T'^ (we denote T = M/27rZ), for some d Cz N. The latter can be motivated by 
the fact that numerical simulations of (1.6) are mainly based on pseudo-spectral 
schemes and thus naturally posed on T"*, see e.g. [2, 3]. We then consider the initial 
value problem for the slightly more general NLS 

(1.6) iedtu' + —Au^ ^ Xelu'l^^u' ; u%0,x) = ul{x), 

where cr g N*. Although we obtain the most precise results (concerning the geom- 
etry of resonances, in particular) in the case of the cubic nonlinearity (cr = 1), we 
are in fact able to rigorously justify WKB asymptotics also for higher order nonlin- 
earities. We assume that (1.6) is subject to an initial data Uq, which is assumed to 
be close (in a sense to be made precise in §6) to superposition of highly oscillatory 
plane waves, i.e. 

(1.7) «g(x)« 5]a,(x)e-^■-/^ 

J&Jo 

where Jq C Z is a (not necessarily finite) given index set. In the Euclidean case 
we allow for wave vectors kj S K'', whereas on = T"* we impose kj G . 
Moreover, in the latter case, we choose ttj to be independent of x e T'', so that (1.7) 
corresponds to an expansion in terms of Fourier series (with S N). The case 
of x-dependent a^'s on T'' could be considered as well, by reproducing the analysis 
on M"^. We choose not to do so here, since it brings no real new information. 

In particular, for a; € T (the one-dimensional torus), our analysis leads to a 
remarkably simple approximation. 

Theorem 1.1. For x G T, consider (1.6) with a = 1. Suppose that the initial data 
are of the form (1.7) with kj = j G Z and {Q:j)j G £^{7j). 

Then for all T > 0, there exist C = C{T) and Eq > 0, such that for all e g]0, Eq], 
with l/e gW , it holds 

sup \\u^{t) ~ ul {t)\\ Ce, 

where the approximate solution it|pp is given by 

<ppit,x) = ^a,e''^^*^''''-l"^■l')e*^^"-4^■'*^/^ and M = ^ |afep. 
jez, fcez 
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We see that at leading order, the nonUncar interaction shows up through an 
expHcit modulation at scale 0{1). It is well known that the one-dimensional cubic 
Schrodinger equation is completely integrable (see [LS, 24] for the periodic case). 
However, this aspect does not play any role in the proof of Theorem 1.1, which in 
itself can be seen as a consequence of the more general result stated in Theorem 6.5. 
On the other hand, several aspects in the discussion on possible phase resonances 
and the creation of amplitudes seem to be specific to both properties d = 1 and 
a = 1 (see §2 and §3). 

In order to prove Theorem 6.5, and henceforth also Theorem 1.1, we need to 
set up a rigorous multiphase WKB approximation for solutions to (1.6). To this 
end, there are essentially two steps needed in our analysis. First, we detail the 
approach sketched above by examining the possible resonances between the phases, 
and analyzing the evolution and/or the creation of the corresponding profiles Uj. 
The second step then consists in making this approach rigorous: we construct the 
profiles Oj, and show that the obtained ansatz is a satisfactory approximation of 
the exact solution u'^, up to 0{e) in a space contained in L°°{A4). As it is standard, 
we prove in fact a stronger stability result: Starting from any approximate solution 
u|pp constructed on profiles, we show that, for any initial data close (as e goes to 
zero) to w|pp|^^Q, there exists an exact solution which is close to u|pp, on some time 
interval independent of e (which, for e small enough, may be chosen as any finite 
time up to which M|pp is defined). 

In the case of a single oscillation only, it suffices to multiply u'^ by e"*"^/^ to 
filter out rapid oscillations, see [G]. In the case where several phases are present, 
this strategy obviously fails. To overcome this issue, a fairly general mathematical 
approach, which has proved efficient in several contexts (see e.g. [15, 27, 2G]), 
consists in working in rescaled Sobolev spaces, usually denoted by H^, for s > 0. 
These are the usual Sobolev spaces, where derivatives are scaled by e, in order to 
account for the spatial oscillations at scale e. More precisely, if s € N, 

\\f\\h--= E wi^drnh- 

\a\^s 

However, due to the negative power of e in the associated Gagliardo-Nirenberg in- 
equalities, this technique usually demands to construct approximate solution with 
a high order of precision (see [14] for a closely related study on the interaction 
of high-frequency waves in periodic potentials). Another, more sophisticated, ap- 
proach consists in filtering out the rapid oscillations in terms of the free evolution 
group, as in [29]. In the present work though, we shall use a simpler approach, which 
allows us to justify the multiphase weakly nonlinear WKB analysis in a remark- 
ably straightforward way. This approach relies on the use of Wiener algebras, as 
introduced in [20], and further developed in [22, 4, 9]. This analytical framework is 
particularly convenient in the case of plane waves, but could probably be extended 
to more general situations, up to some geometric constraints on the phases. How- 
ever, the first step of the analysis, i.e. describing all possible resonances, becomes 
much more intricate, see e.g. [21, 19]. 

As well shall see during the course of the proof, the use of Wiener algebras has 
several advantages on the technical level. We point out that this framework makes 
it possible to justify the WKB approximation with an error estimate of order C(e) 
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without constructing correctors (which would have to be of order e or even smaller, 
when working in spaces, see e.g. [15, 27], or [7, 14] in the NLS case). 

1.4. An application to instability. As an application of the semi-classical anal- 
ysis for (1.6), we recover the main result in [cS] (see also [-5]), concerning NLS in 
the periodic case. This result has been established in the case d = 1, and is hereby 
extend to higher dimensions. We also propose a variation on a result in [2-5] (see 
assertion 3 in the theorem below). 

Theorem 1.2. Let d ^ 1, cr e N* and X e {±1}. Fix s <0. 

1. For all p > 0, we can find a solution u to 



for some constant c > independent of p and 5. In particular, the solution map 
fails to he continuous as a map from H^iT'^) to H''(T'^), no matter how close to 
— oo the exponent k may be. 

2. Suppose (7^2. For any p > and S > there exist smooth solutions u, u of 
(1.8) such that u(0) — u{Q) is equal to a constant of magnitude at most 5, and 



for some constant c > independent of p and S. 

3. For any t 0, the flow-map associated with (1.8) is discontinuous as a map from 
L'^{T^), equipped with its weak topology, into the space of distributions (C°°(T'^))* 
at any constant ao G C \ {0} C L'^iT'^). 

We show in §7 that the above instability result can be viewed as a consequence 
of multiphase weakly nonlinear geometric optics. The first two assertions are an 
extension of the results in [8] , so we shall not comment on their meaning, and refer to 
the discussion in [8]. We invite the reader to consult [25] for a stronger instability 
result in the one-dimensional case: indeed, when d = a = 1, the author shows 
the third point in the above statement for any ao e i^(T) \ {0}, not necessarily 
constant. 

1.5. Structure of the paper. We first study in detail the case of the cubic nonlin- 
earity (cr = 1). In §2, we consider the set of resonant phases, and in §3, we analyze 
the corresponding amplitudes. The case of higher order nonlinearities is treated 
in §4. In §5, we set up the analytical framework, with which a general stability 
result (of which Theorem 1.1 is a straightforward consequence) is established in §6. 
Theorem 1.2 is proved in §7. Finally, in an appendix, we sketch how the previous 
semi-classical analysis can be adapted to more general sets of initial plane waves 
(including generic finite sets of wave vectors). 

Acknowledgments. The first author wishes to tliank Thierry Colin and David Lannes 
for preliminary discussions on this subject. 




with \\u{0)\\h=('ii1) < P, such that for all S > 0, there exists u solution to (1.8) with 

||u(0) -m(0)||h.(t<') < S, 
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2. Analysis of possible resonances in the cubic case 

In this section, we show that when tr = 1, the set of relevant phases can be 
described in a fairly detailed way. 

2.1. General considerations. We seek an approximation of the form 

je.J 

where here and in the following J C Z denotes the index set of relevant phases 
4>j (yet to be determined). Note that using J is only a renumbering, so that 
j ^ k ^ (j)k- In the case x gT^, one simply drops the dependence of aj upon 
X. In general Jq C J, i.e. we usually need to take into account more phases in (2.1) 
than we are given initially. 

As a first step we need to determine the characteristic phases (j>j {t, x) G K. For 
plane- wave initial data of the form (1.7) we are led to the following initial value 
problem 

dt<t>o + ^ I V0jf = ; <t>j (0, x) = ■ X, 
the solution of which is explicitly given by 

(2.1) (/)j(<,a;) = Kj-a;-||Kjf. 

Recall that for x gM.'^, we assume Kj £ R'^, whereas in the case i £ T'^, we restrict 
ourselves to Kj S Z''. Of course, these phases <j>j remain smooth for all time, i.e. 
no caustic appears. 

In the cubic case cr = 1, the set of resonances leading to the phase (f>j is therefore 
given by 

Ij ^ {{k,£,m) e ; Kk ~ Kg + k,^ ^ Kj, l^fc^ - l^fP + l^mP = l^jP}, 
and the corresponding resonance condition (1.5) becomes 

(2.2) (k£ - Kjn) ■ (Ki - Kk) = 0. 

As we shall see, this condition provides several insights on the structure of four-wave 
resonances. 

2.2. The one-dimensional case. For d = 1 the condition (2.2) implies that if 
{k,i,m) S /j, then ki ~ Kj„, or ki — Kk- Therefore, when d = 1, the set Ij is fully 
described by: 

I, ^ {{,],£,£), i£,£,j) ; J}, 
no new phase can be generated by a cubic interaction. 

In higher dimensions, however, the situation is much more complicated and 
heavily depends on the number of initial modes. 

2.3. Multi-dimensional case d ^ 2. We start with the simplest multiphase situ- 
ation and proceed from there to more complicated cases. Eventually we shall arrive 
at a geometric interpretation for the generic case. 
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2.3.1. One or two initial modes. If we start from only two initial modes, (J^/o = 2, 
the resonance condition (2.2) implies that the cubic interaction between these two 
phases cannot create a new characteristic phase. In other words, exhibits at 
most two rapid oscillations at leading order. Recalling that = is an admissible 
phase, the case of a single initial phase jJJo = 1, is therefore included (if one of the 
initial amplitudes is set equal to zero). We want to emphasize that the case of at 
most two initial phases is rather particular, since (2.2) implies that the situation is 
the same for all spatial dimensions d ^ I. 

Remark 2.1. In addition, the fact that two phases cannot create a new one extends 
also to higher order (gauge invariant) nonhnearities f{z) = Ajzp'^z, for ct G N, 
0-^2, see §4. 

2.3.2. Three or four initial modes. This case can be fully understood by the follow- 
ing geometric insight, already noticed in [10]: 

Lemma 2.2. Let d ^ 2, and k,£,m belong to J. Then, (nk, Kg, k„i) G Ij pre- 
cisely when the endpoints of the vectors k^, Kg, Km, Kj form four corners of a non- 
degenerate rectangle with K( and Kj opposing each other, or when this quadruplet 
corresponds to one of the two following degenerate cases: (^A; — , — ^i) ; 

i^^k — — ^j) ■ 

Remark 2.3. In the degenerate cases, no new phase is created. 

Proof. We recall the argument given in [10], by first noting that the relations be- 
tween (Kj, Kk, Kg, Km) formulated in (2.1), are equivalently fulfilled by (k^ — k, ki — 
K,Km — K,Kj — k), for any k G M.'^ (resp. k E Z'^). This is easily seen by ex- 
panding the second relation in (2.1) and inserting the first one. Thus, choosing 
K = Kj, it therefore suffices to prove this geometric interpretation for Kj — 0, which 
consequently shows: Kk + = such that Kk ■ k„i = 0, by the law of cosines. □ 

In summary, we conclude that three initial (plane-wave) phases create at most 
one new phase, such that the corresponding four wave vectors form a rectangle. 
When the initial wave vectors {tjljeJo are chosen such that their endpoints form 
the four corners of a rectangle, no new phase can be created by the cubic nonlin- 
earity and u'^ exhibits only four rapid oscillations. We close this subsection with 
two illustrative examples. 

Example 2.4. Let d = 2. Consider ki = (0,1), K2 = (1,1) and K3 = (1,0). The 
cubic interaction creates the zero mode (p/^ = 0. 

Example 2.5. Again let d = 2, with now ki = (1, 1), K2 = (1, 2) and K3 = (3, 2). In 
this case, we create a non-zero phase 04, with corresponding wave vector K4 = (3, 1). 

The geometric insight gained above then directly leads us to the following de- 
scription of the resonant set Ij in the general case. 

2.3.3. The general case. We are given a countable (possibly finite) number of initial 
phases {4>j}jeJQ with corresponding wave vectors {Kjjjgjg. From the discussion of 
the previous paragraph it is clear that there are two possible situations: 

(a) Either, it is impossible to create a new rectangle from any possible subset 
Jo G Jo, such that [JJo = 3. If so, then no new phase can be created. This 
is the generic case. 
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(b) Or, starting from an initial (finite or countable) set 5*0 = {'«j}je,/oj we may 
obtain a first generation Si = {^jljeJi with Jq C Ji (i.e. 5*0 C Si) in 
the following way: we add to 5*0 all points k e M'', such that there exist 
Jo C Jo with ft Jo = 3; and such that {^j} U {k} is a rectangle. Note that, 
if Jo C Z'', then Ji C Z''. By a recursive scheme, we are led to a (finite 
or countable) set S which is stable under the completion of right-angled 
triangles formed of points from this set, into rectangles. Furthermore, if 
So C Z*^, then S C Z'K 

Example 2.6. As already seen, the simplest examples for possibility (a) are the cases 
tt Jo ^ 3 when the triangle formed by the endpoints of the considered wave vectors 
has no right angle, or (jJo = 4, where the four initial phases are chosen such that 
their corresponding wave vectors {'tjljeJo already form the corners of a rectangle. 

From a finite number of initial phases, possibility (b) may lead to a finite as well 
as to an infinite set J. Even for d = 2, we have: 

Example 2.7. In the plane R^, start with 

Jo = {(-1,1), (0,1), (0,0), (1,0)}. 

The first generation is then 

Ji = {(-1, 1), (0, 1), (1, 1), (-1, 0), (0, 0), (1, 0)} = Jo U {(1, 1), (-1,0)}, 

and the second one is 

J2 = JiU {(0,2), (0,-1)}. 
One easily sees that this generates J = Z^. 

As a conclusion, the set of phases {(/)j}jgj may be finite or infinite, but has the 
following property. 

Proposition 2.8. Let (7 = 1, and consider any triplet of wave vectors from S = 
{Kj}jizj. Then, either the corresponding triangle has no right angle, or the fourth 
corner of the associated rectangle belongs to S . 

3. Analysis of the amplitude system in the cubic case 

From the previous section, in general we have to expect the generation of new 
phases by the four-wave resonance. However, it may happen that not all of them are 
actually present in our approximation (2.1), since the corresponding profile aj(t,x) 
has to be non-trivial. 

Indeed, if we plug the ansatz (2.1) into (1.6) the terms of order 0(1) are iden- 
tically zero since all the (/)j's are characteristic. For the 0{e) term, we project on 
the oscillations associated to (pj, which yields the following system of transport 
equations: 

(3.1) Vj e J, dtOj + Kj ■ Woj = —iX ^ akaiOm ; aj{0,x) = aj{x), 

{k,e,m)eij 

with obviously Voj = in the case where x E T'^. In the following we will perform 
a qualitative analysis of the system (3.1), postponing the rigorous existence and 
uniqueness analysis to §5.4. Having in mind the discussion from §2 we distinguish 
the case d = 1 from the case d ^ 2. 
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3.1. The case d = I. Let j e J, and recall that Ij is particularly simple in d = 1: 
Using this, (3.1) simplifies to 

(3.2) {dt + Kjdx) aj = —'2iX''^^\a(faj + i\\ajfaj ; aj{Q,x) = aj{x). 

In particular, the evolution of a zero profile aj = is necessarily trivial, that is 
aj{t,x) = 0. This non-generation of profiles leads to the same conclusion as §2.2: 
No new mode can be created, if it is not present initially (and the reason is the 
same as in §2.2: Uj factors out in (3.2) just because for any (£1,^2,^3) S Ij, we 
have ii = j or ^3 = j). We shall see that the multi-dimensional situation is quite 
different but first examine the situation for .t S T and a: € M in more detail. 



3.1.1. The case a; G T. In this case, we readily obtain that |ajp does not depend on 
time. This is due to the fact that (3.2) yields: idtaj G M.aj and hence 9t|ajP — 0, 
for all j G Z. In particular we get that 

The conserved quantity M corresponds to the total mass of the exact solution u'^. 
Using this, we rewrite (3.2) as 

^flj = -iX (2M- 
which yields an explicit formula for the (global in time) solution 

We observe that in the case of the one-dimensional torus, the interaction of the 
profiles aj is particularly simple. Nonlinear effects lead to phase-modulations only. 

3.1.2. The case cc € M. Here, in contrast to the situation on T, the modulus of aj 
is no longer conserved, since we can only conclude from (3.2) that 

{dt + Kjd,)\aj\^ = 0, 

and thus 

\aj{t,x)\^ = \a.j{x-tKj)\^. 
In particular we readily see that for all j G J we have 

(3.3) \\aj{t)\\L2 = \\a,\\L2, Vi G R. 

Moreover, we still have an explicit representation for the solution of (3.2) in the 
form 

(3.4) aj{t,x) = aj{x - tKj)e^'^'^*-''\ 

for some real- valued phase 5"^, yet to be computed. In view of the identity 
{dt + Kjd^) aj{t, x) = iaj{x - tKj)e*^'^J^*'^^ {dt + Kjd^) S.j{t, x), 
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equation (3.2) implies 

{{dt + Kjdx) Sj{t, x)) aj{x — tKj) =\ j —2 |a^(a; — tK^)P + \aj(x — tKj)\'^ 

X aj(x — tnj). 
One easily sees that it is sufBcient to impose 

dt {Sj{t,x + tKj)) = -2X^\ae{x + t{Kj - k^))]^ + X\aj{x)\'^, 

which yields 




^3^^ Sj{t,x)^-2\l ( \aeix + iT-t)K,^TK,))\'dT 

— tX\aj{x — tnj)]'^ . 

This formula, together with (3.4) describes the modulation of the profile aj{t,x). 
As in the case of the torus, amplitudes are transported linearly. Only the (slow) 
phases Sj undergo nonlinear effects, which are more complicated as before but still 
explicitly described in terms of the initial data. 

3.2. The case of one or two modes for d ^ 1. We have already seen in §2.3.1 
that the case of two initial modes is special, since we get a closed system for all 
d ^ 1. Indeed if we start from two phases and two associated profiles, say Uj and 
ae, the system (3.1) simplifies to: 

dtttj + Kj ■ \7aj = —iX (|ajp + 2|afp) Uj, ; aj{0, x) — aj{x), 

dttti + Ke -Vae = -iX (2\aj\^ + \a£\'^) Qi, ; ai{0, x) = a£{x). 

Note that if initially one of the two profiles is identically zero, it remains zero for 
all times and hence, we are back in the situation of a usual single-phase WKB 
approximation. In particular we compute explicitly for: 

• Two modes, on T'': 

• Two modes, on M**: 

aj{t,x) = ajix ~ tKj)e-'^(^^« |a,(x+(r-t)«,-rK,)| = rfr+t|a,(x-t«:,)|^)^ 
ai{t,x) = ai{x - tHi)e-'^(^^o |a,(x+(r-t)«f-r«,)| = rfr+t|af(a.-t«,)|^) ^ 

Again, these solutions exhibit (nonlinear) self-modulation of phases only, and exist 
for all times t e R, a property which is a-priori not clear in the general case. 

3.3. Creation of new modes when d ^ 2. A basic difference between the one- 
dimensional case and the multidimensional situation is that the conservation law 
(3.3) does not remain valid when d ^ 2. However, we are still able to prove that 
the total mass is conserved. 

Lemma 3.1. For any solution to (3.1) it holds 

(3.6) j^J2\Mmh = o. 

je.J 
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Proof. The assertion follows from the more general identity 

^ (a* + • V) |a,f = 0, 
jeJ 

since, by definition we have 

{dt + Kj ■ V) Iflj p = Im A ^ ^ ajUkaeam ■ 
jeJ \ jeJ {k.e.m)eij / 

This sum is zero by symmetry, since for each quadruplet (j, k, i, m) € J* of indices 
the quadruplet (j, m, £, k) is also present, as well as the other six obtained by 
circular permutation (at least in the nondegenerate case mentioned in Lemma 2.2; 
adaptation to the degenerate case is obvious). These are the only occurrences of 
the corresponding rectangle of wave numbers, and they produce the sum 

2 (ajOfeaf a„ + akaeamOj + aeamajak + amajOkai) = 8 Re {ajakaiam) , 

which is real. We consequently infer 

dt^\aj{t,x + tKj)\'^ = 0, 

and thus also 

□ 

Let us now turn to the possibility of creating new profiles by nonlinear inter- 
actions (note however that the conservation law (3.6) gives a global constraint on 
this process). To simplify the presentation, we assume d ~ 2. The creation of 
new oscillations in the general case d ^ 2 then follows by completing elements in 
with (0, . . . , 0) G K''"^ and analogously for the situation on T*^. Consider the 
geometry associated to Example 2.4: We thus have (on T'' or R'^) 

idtoo = A E akaettm- 
Recall that (fc,^, m) G /o if and only if 

Kfe-K£ + Km=0 ; - + |k„iP = 0, 

which obviously implies • Km = 0. Such a possibility occurs in two cases: 

• K/j = or Km — 0. 

• iKk,Km) = (ki, K3) or (Kfc, Km) = (^3, Ki) and hence k/ = K2. 
From these various cases, we infer 

idtao = A (|aop + 2|aip + |a2p + 2|a3p) ao + 201^203. 

Consider three non-vanishing initial oscillations, such that aia2a3|t=o 7^ 0. Thus, 
even if ao|t=o = 0, we have 9tao|t=o 7^ 0; ^-nd this (non-oscillating) fourth mode is 
instantaneously non- vanishing. 
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4. Higher order nonlinearities 

4.1. Analysis of possible resonances. So far wc were only concerned with four- 
wave interactions corresponding to cubic nonlinearities, i.e. a = 1 in (1.6). In 
general though, the set of resonances associated to a (gauge invariant) nonlinearity 
of the form /(z) = Ajzp'^z, cr G N, are defined by 

k=l k=l 

As in Section 2, the set of wave vectors {Kj}ji=j is constructed by induction, starting 
from an a finite or countable set {Kj}j^jg, to which we first add a vector k when 
there exist Kg-^, . . . , k^^^^^ € Jo such that 

(4.1) E (-1)''+V£j' = 

k=l 

we then set k = J^k^^^ {~^)'^~^^ . The same iterative procedure as in §2.3.3 leads 
to the following analogue to Proposition 2.8: 

Proposition 4.1. Let a ^ 2, and consider any {2a + l)-tuple (k^^, . . . , k^2<t+i) 
wave vectors from S = Then, either the relation (4.1) is not satisfied, or 

the vector kj ~ '^'k^^ i^^)'^^^ i^ik belongs to S. 

Remark 4.2. It is worth noting that, even if we only have very poor information 
on the set of wave vectors {Hij}j^j, it is however a subset of the group generated 
by the initial set {Hij}j^jg. 

The profile equations, analogue to (3.1), are then, for all j E J: 

(4.2) dtOj + Kj ■ Vuj = —iX a^iftfa • ■ ■'^t2a+i ! aj{0,x) ~ aj(x). 

(£i,...,«2^+i)e/j 

4.2. The case of two modes. Similar to the situation for a = 1, the case of 
only two initial modes is rather special. Indeed, the fact that two phases cannot 
create a new one extends also to higher order nonlinearities. In order to explain the 
argument, consider first a quintic nonlinearity, corresponding to cr = 2. To obtain 
a nonlinear resonance, the wave vectors need to satisfy 

l^k K/i -\- K„i Kp -\- Kq = Kj, 

where k,£,m,p,q € {ji, ^2}, Ji, ^2 £ J- First, if ji (or J2) appears at least twice on 
the left hand side, with at least one plus and one minus, then the cancellation re- 
duces the discussion to the one we had about the cubic nonlinearity. Hence, no new 
resonant phase can be created in this case. The complementary case corresponds, 
up to exchanging ji and j2 , to 

The above relations yield 

3Kjj — 2Kj2 = Kj ; 3|kjJ^ — 2|Kj2p = 
Squaring the first identity and comparing with the second one, we infer 

6|«:ji - 1^32? = 0- 



2cr+l 

fe=i 



\k+l 
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Therefore, no new resonant phase can be created by the quintic interaction of two 
initial resonant plane waves. 

Consider now the general case where cr > 2: The same argument as above shows 
that the only new case is the one where all the plus signs correspond to one phase, 
and all the minus signs to the other: 

(cr + - CTKj, = ; {a + l)\Kj,\^ -a\K,jJ^ = \nj\^. 

Squaring the first identity and comparing with the second one, we infer 

a{a + l)\Kj^ - = 0. 

We conclude as above, and obtain the following result: 

Proposition 4.3. Let a e N*, and let ki,K2 € be such that ki ^ H2- To these 
wave vectors are associated the characteristic phases 

(/)j(t,a;) = • a; - j==l,2. 

Then, these two phases can not create no new phase by {2a+l)th-order interaction: 
the set 

2(7+1 

k=l 

2(7+1 



k=l 



is reduced to {ki,K2}- 

In view of Proposition 4.3, the system (4.2) becomes a system of two equations, 
which can be integrated explicitly, as in [S, Remark 3.1]: 



dfUj + Kj ■ Voj = —iX 



n=0 



a + 1 \ f a 
n I \ n 



n \ n 



(4.3) ^ . . . 
dtai + ti( ■ Vflf = — iA ( 

)i=0 ^ 

In the case of T^, wc find for instance 

/ 

aj{t) = aj exp j — iA< 

(4.4) ^ 

aeit) = exp j —iXt 



W (7 
(7+ 1 \ / CT 

71 I \ n 



\2(7 — 2n\ \2n 



\2a~2n\^ \2n ^ 



\2a~2n\„ \2l 



|„ |2(T — 2n|„ |2n 



\ ri=0 

In the case of M'', the formula is more intricate and we shall omit it. 

Apart from the two- phase situation, the results for of Section 2.3.3 on resonances 
do not carry over to the general case tr ^ 2 in any straightforward manner. Even in 
space dimension d = 1, the resonant sets cease to be as simple for cr ^ 2, provided 
that one starts with at least three modes. 
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Example 4.4. Consider the quintic case ct = 2 in d = 2 spatial dimensions. As we 
have seen above a resonance for such a quintic nonhnearity appears if and only if 

l^k K£ -\- K„i Kp -\- Kq = Kj, 

We can pick for instance three initial phases of the form 

Ki = (-1,0) ; «:2 = (0,0) ; Kg = (2,0). 

For k = 1, i = p = 2, m = q = 3, we have a resonance, creating K4 = (3, 0), whereas 
in the case ct = 1, no resonance occurs between the phases with wave vectors ki, 
K2 and K3. This example shows that the geometric characterization of four- wave 
resonances given in §2.3.2 does not export to the case of six- wave resonances: ki, 
K2, K3 and K4 all belong to the line X2 = 0- 

Example 4.5. Consider the same example as above in d = 1. i.e. pick three initial 
phases of the form 

Ki = — 1 ; K2 = ; K3 = 2 

and create a resonance K4 — 3 ioi k = 1, i = p = 2, m = q = 3. This is in sharp 
contrast to the case cr = 1, where no new phases can be created in d = 1. Moreover, 
a non- vanishing amplitude 04 is effectively generated: 

{dt + «49.)a4 = - 3a (|ai|4 + jagl'* + lasi' + 4(|aina2|' + \a2\^\a3\^ + Wsl^Wil^)) 

— i6Xaia2^al — 6zAai|a2pa3. 
We see that we may have a4|t^o = 0, but 

(9ta4)|j^o = (-iBAaioi'^ag - 6iAai|a2pa3) i^^^ ^ 0, 
showing the appearance of a non-trivial 04 for t > 0. 

Despite this lack of knowledge concerning the precise structure of possible res- 
onances for higher order nonlinearities, we shall sec that we arc able to prove the 
validity of WKB approximation even in this case. 

5. Analytical framework 

We now present the analytical framework needed for the rigorous justification of 
a multiphase WKB approximation. 

5.1. Wiener algebras. On M = T^, we consider the usual Wiener algebra of 
functions with absolutely summablc Fourier series: 

Definition 5.1 (Wiener algebra on A4 = T^'). Functions of the form 

f{y) - bkeJ^"""^ with Kk e Z"- and bu G C, 
fcez 

belong to WiT'^) if and only if {bk)kez S i^C^). We denote 

\\f\\w^Y.\bk\. 

In the sequel, when x G T'', we consider initial data for (1.6) which are of the 
form f{x/e), with / G WiT'^) and G N*. 
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Lemma 5.2. Let f belong to WiT"^). Then, for all e > such that e ^ G N*, we 
have /(•/£) e WiT'^), and 

\\fi-/e)\\w^\\f\\w. 

For = K."^, the framework is a bit different. Define the Fourier transform by 

With this normalization, we have f{x) ~ Tf{—x). Following ['2i)] and [')], we 
use on W'- two different Wiener-type algebras: For the exact solution we use VF(]R''), 
i.e. the space of functions with Fourier transform in L^{W^). and for the profiles, 
we use A(E^), the space of almost periodic VK(K'') -valued functions on R'', with 
absolutely summable Fourier series. We also set A(T'^) = P^(T''), equipped with 
the same norm. 

Definition 5.3 (Wiener algebra on = M''). We define 

WiW") = {/ e 5'(M'^;C), WfWw WIWlhr^) < ^} ■ 

Functions of the form 

f{x,y) ^Y^k{x)e""'-y, with e and bk <E W{R'^), 
fcez 

belong to A(]R'^) if and only if 

||./||a := Y ll^'^-ll^' = IMlhm") < oo. 

In the sequel, when x G R'^, we consider initial data for (1.6) which are of the 
form f{x,x/e), with / e A(R''). Again, we have 

Lemma 5.4. Let f e A(R'') and e > 0. Then /(•, -/e) G W{R'') and 

\\f{-r/e)\\w^\\f\\A. 
Proof. We simply have, when f{x,y) = J2kez^k{x)e^'^''"'': 

\\fi-r/£)\\w = II Xl^'^('~^'=/^)ILi(R'') ^ ^ \\hi-~Hk/e)\\LHM'i} = Y W^kh^R")- 
fcez feez feez 

The last term is, by definition, ||/||a- D 

Denote (in the periodic setting as well as in the Euclidean case) 

U'it) = e'^^A_ 

The following properties will be useful (see [9], and also [20, 22, 4]). 
Lemma 5.5. Let M^T"^ or M''. 

1. W{A4) is a Banach space, continuously embedded into L°°{A4). 

2. W{M) is an algebra, in the sense that the mapping {f,g) fg is continuous 
from W{A4)^ to W{Ai), and moreover 

yf,geW{M), WfgWw^WfWwMw. 

3. Lf F : C ^ <C maps u to a finite sum of terms of the form uPyfl, p,g Cz N, then it 
extends to a map from W{A4) to itself which is uniformly Lipschitzean on bounded 
sets ofW{M). 

4. For all t G M, U^{t) is unitary on W{M). 
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5.2. Action of the free Schrodinger group on W(Ai). As it is standard for 
solutions to the equation 

iedtw^ + —Aw^ — F^, 
we will consider the corresponding Duhamcl's formula 

Jo 

In view of this representation formula we first need to study the action of the 
free Schrodinger group U^{t) on W{A4). 

5.2.1. The case Ai = T^. The action of U^{t) on Fourier series on T'* is well 
understood. For Efeez^fee"""'^ S W{T''): 



(5.1) 



rt/(2e) 



In view of Duhamcl's formula, we will use the following 

Lemma 5.6. LetT >0, uj e Z, k e , and b, dtb e L°°{[0, T]). Denote 

1. We have e C{[0,T] x T'^) and 

II-D1Il-([o,t]xT'') / \b{t)\dt. 

Jo 

2. Assume u) ^ j/tp. Then there exists C independent of k, uj and b such that 

Ce 

ll-0^|lLoo([o,T]xT'i) Tj— 1^ 1 {\\b\\L^{[0,T]) + II ^^t^ll L°= ([0,T])) ■ 

Proof. In view of the identity (5.1), we have 

D'{t,x) ^ f\(r)e'^-^^^-'^-/i^^)e-'\^\'^t-^y(^^)dT 
Jo 



The first point is straightforward. Integration by parts yields, since by assumption 
|k|2 - w e Z \ {0}: with (j){t, x) ^ k ■ x - \K\^t/2, 







V l^l"' — UJ 

\k.\--ujJq / 
The lemma then follows easily. □ 
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5.2.2. The case M. = . The Euclidean counterpart of Lemma 5.6 is a little bit 
more delicate: 

Lemma 5.7. Let T > Q, uj k <eW^ , and b <E L°°{[0,T];W{M.'^)). Denote 
D'{t,x) := U'it - t) (6(T,x)e*'"-^/^-''^^/(2c)^ 

1. We have e C{[0,T];W{R'')) and 

\\D'\\L^ao,T];W) ^ I \\b{t,-)\\^dt. 





issume ljJ ^ and dtb, Ab £ L°°([0,T]; W). Then we have the control 
Ce 



\\D''\\l°-(10,T];W) ^ ||^|2 _ ^1 (II^IIl~([0,T];IV) + II^^IIl~([0,T];H') + ilO. T]-W) ) ' 

where C is independent of k, uj and b. 
Proof. By the definition of U'^{t), we have 

D'{t, e) = e-'^(*-^)l«l'/2 6 (r, ^ - ^) e~'^^/^^''>dT. 

Setting rj = — n/e, we have 





where we have denoted 

9 = e f] 



K UJ , ,9 „ IkP — UJ 



e - ^ ' e 



The first point of the lemma is straightforward. To prove the second point, integrate 
by parts, by first integrating e''^®^/^: 

D%t, - -^/^'^'bir, ry) [ + I e"'l^ (^z^fo (r, ^) + d^b (r, r;)) dr. 

We infer, if 6, (9t6, Afe e L°°([0, T]; W^): 

sup ||i5''(t)||ii < (||6||l-([0,T];L1) + l|A6||i=o([o,T];Li) + \\dtb\\L--{[a.T]:L^)) ■ 

te[0,T] |P2| ^ ■ ■ / 

This yields the second point of the lemma. □ 

5.3. Construction of the exact solution. As a preliminary step in establishing 
a WKB approximation we first need to know that (1.6) is well posed on W{A4). 

Proposition 5.8. Consider the initial value problem 

(5.2) iedtu" + —Au^ ^ \e\u^\^''u^ ; u''(0, x) = u^(a-), 

where a £ W, A e M, and x € M, with either M = K'^, or M = T'', in which 
case G N*. // Wq € W{M.), then there exists > Q and a unique solution 
e C{[Q,T^];W{M)) to (5.2). 
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Remark 5.9. At this stage, the dependence of upon e is unknown. In particular, 
might go to zero as e — > 0. The proof below actually shows that if Uq is uniformly 
bounded in W{A4) for e g]0, 1], then > can be chosen independent of e. This 
case includes initial data (1.7) which we consider for the WKB analysis. 

Proof. Duhamel's formulation of (5.2) reads 

u%t) = U'{t)uf, -i\ f U'{t - r) {lu'l^^u^iT)) dr. 
Jo 

Denote by $^(u^)(t) the right hand side in the above formula. From Lemmae 5.5, 
5.6 and 5.7, we have: 

||<i>-(u^)(t)||^^ ^ IKIIw + |A| r \\u^T)f^+'dr, 

Jo 

and if ||w^||l^([o,t];W)i ll^i^l|L°^([o,T];H') ^ then there exists C = C{R) such that 

\\^%u^m-^^{u^){t)\\^^C{R) f\\u^{T)~u%r)\\^.dT, Vte[0,r]. 

Jo 

A fixed point argument in 

\ueC{[(),T]-W{M)), sup \\u{t)\\w ^2\\ul\\w\ 

{ tG[0,T] J 

for T ^ > Q sufficiently small then yields Proposition 5.8. □ 

5.4. Construction of the profiles. In order to justify our multiphase WKB anal- 
ysis, we first need to establish an existence theory for the system of profile equations. 
To this end, for all cr € N*, we rewrite the system (4.2) in its integral form: 

(5.3) V j € J, aj{t,x) = aj{0,x ~ tKj) ~ iX / A/'a(a, . . . , a)(r, a; + (r — t)Kj)c?T, 

Jo 

where, for a^^^ = (a^^')jgj, . . . , a^^'^+^' — (aj^'^^^^)jgj, we define the nonhnear 
term Afa by: 



VjeJ, AA.(aW,...,a(2'^+i))= ^ 



(1)_(2) _(2o-) (2cr+l) 

a^^ ... a^^^ a^2<T+i 



It is clearly linear with respect to its arguments with odd exponents, and anti-linear 
with respect to the others. We prove in Lemma 5.11 below that it is in fact well 
defined and continuous on E{M), for M =T''- or M ^M'^: 

Definition 5.10. Define 

E{R'') = {a= {a,),ej \ {^,)jeJ € t\j ; L^W'))} , 

equipped with the norm 

\\a\\E(S.<^) = X! Il^jllii- 
Set also E{T'^) = £^{J), equipped with the usual norm 
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Note that E simply represents, via an isometric correspondence, the family of 
coefficients of functions in A (up to the choice of the wave numbers kj in the case 
of R'^): 

f{x,y) ^aj(a;)e*'"^'^ G A(R'') iff a e E{R'^), 
and then ||a||£; = ||/||a- The same holds for M ~ T''. 

Lemma 5.11. Let ct G N*. For Ai = M'* or Ai = T'', the nonlinear expres- 
sion Na defines a continuous mapping from E{AiY'^^^ to E(Ai), and for all 
a(i),...,a(2<T+i) g E{M) 



) ,,a(2-+i) 



Proof. We consider the case Ai = M."^, since = T'' is even simpler. In order to 
bound 



E 

je.J 



(l....,l2a+l)elj 



E 

je.i (<'i,.../2<T+i)e/j 



(2: 



we use Young's inequality and observe that, once j, £1, . . . , £20- ^ire chosen, £2^+1 

is determined (since K£2<t+i = '^j ^ 12kLii^^)'^^^ '^^k^ ^^'^ n ^ m ^ k„ 7^ Km), so 
that 



AA^(a(i',...,a(2-+i)) 



E 



i)eJ2-+i 

















Ll 



which gives the desired result. 



□ 



This consequently yields the following existence result for (5.3), where here and 
in the following we denote {nf' = 1 + 

Proposition 5.12. Let ct G N*, and Ai=m.'^ or A4=T'^. 

For all a = {aj)j^j € E(A4), there exist T > and a unique solution 

t^a{t) - {a,{t)),e.J e C{[0,T],E{Ai)) 

to the system (5.3), with a(0) = a. Moreover, the following properties hold: 

1. If{{K.j)'a,)j<.j e E{A4.) for some s e N, then {{k,Y aj),^j G C{[0,T], E{M)). 

2. On M ^ R'^, if {{kj^ dl^aj)jeJ e E{R'^), for some (3 e N'' and s e N, then 
{{n,ydSa,),ejeC{[0,T];E{R'i)). 

Proof. The existence result follows from Lemma 5.11 and the standard Cauchy- 
Lipschitz result for ODE's. Concerning the propagation of moments {KjY Oj, we 
again apply a fixed-point argument, estimating nonlinear terms 



(K,rAA.(a«,...,a(2-+i)) 
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as in the proof of Lemma 5.11, via 

^1 + |.^,P = 1+ ^(-1)'=+1|«,J2 
fe=l 

2(7 + 1 2(7 + 1 

E («:^j'<(2t + i) n («4>', 

k=l k=l 

when (^1, . . . , ^2(7+i) G ^j- The last statement of the proposition is concerned with 
the smooth dependence upon the parameter x. This follows by commuting (4.2) 
with dx and using the fact that is an algebra, continuously embedded in 

since then 

-T:\\dxa\\E < \\dxa\\E + C{\\a\\E)\\dxa\\E, 
at 

and a Gronwall argument shows that ||3a;a||£; remains bounded for all t G [0,r]. 
Similarly we conclude for the higher order derivatives, possibly multiplied by weights 
{k,Y- ' □ 

For the particular situation for <t = 1. in c? = 1 and/or the case of only two initial 
phases, we infer a stronger result, thanks to the explicit formulas given in §3.1 and 
§3.2. 

Corollary 5.13. Under the assumption of Proposition 5.12, in the case a = 1, if 
in addition d = 1, then T can he taken arbitrarily large, with aj(t) explicitly given 
by (3.4) and (3.5). Similarly, i/tJJo ^ 2, then T can be taken arbitrarily large. 

Remark 5.14. In the case of higher order nonlinearities, i.e. a ^ 2, Equation (4.3) 
makes it possible to see, via explicit integration (see (4.4) in the case of the torus), 
that if aj,a(, £ W[M), then aj,ai £ C {[Q , oo\,W {M)) . 

6. Rigorous justification of the multiphase WKB analysis 

6.1. Construction of an approximate solution. We start from oscillating ini- 
tial data, given by a profile in A(A^), with = T'' or R'': 

<pp(0,x)^ E".'-(^)^'"^''^'' 

jeJa 

with aj{x) ~ Const, in the case Ai = T''. 

Assumption 6.1. For both 7W = M'' and M ~ T"^ we assume {aj)j^j„ £ E{M). 
For M = we assume in addition 

V|/3| 2, (af aj)je./o e m"), and V|/3| ^ 1, {{k,) a,),ej„ e E{M.''). 

From Proposition 5.12 we know, that these data produce a solution {aj)j^j G 
C{[0,T], E{M)) to the amplitude system and we consequently define the approxi- 
mate solution u|pp by 

(6.1) <pp(t,a;) = E«.'(^'^)e*'*''^^'' 

with (f>j given by (2.1). The sequence {aj)j^j is such that 

{{^3) 5f %-),eJ ^ ^([0, T],E{M)), \P\ ^ 1. 
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We see from equation (5.3) that {dtaj)j^j G C([0, T],E{M)). Wc find (in the sense 
of distributions) 



*£5t<pp + Y^<pp = ^eKppP^Wapp - ^e*-! + e^^2: 



where 
(6.2) 



and the remainder rf takes into account the non-characteristic phases created by 
nonhnear interaction. This means that it is a sum of terms of the form 



where the rapid phase is given by 



2cr+l 



2cr+l 



p=i \ p=i J 



p=i 



and 



2<T+1 



2(7 + 1 

+ E ^-^y^" 
p=i 



In other words. (£i, . . . ,^2ct+i) belongs to the non-resonant set 

N := \ y //. 

With these conventions, we have 

(6.3) rl^ a€,a,,...a,,„^,e^(*^i-'^'=+--'^'-+*Wi)/^ 
(<'i,...,f2^+i)ew 

Estimating in H/^ is straightforward, since {d^aj)j^j € C([0, r],i?) for \j3\ ^ 2: 



(6.4) 



Note that r2 simply vanishes if = T^. In order to estimate rf , we impose the 
following condition on the set of wave numbers 

Assumption 6.2. There exists c > such that for all (^i, . . . , ^2cr+i) G N , 



S ih, . . . 



2(7 + 1 

p=l 



2(7 + 1 

p=i 



^ c. 



Remark 6.3. (i) This assumption is of course satisfied when only finitely many 
phases are created [tJ < c». 

(a) Similarly, this assumption holds for {Kj}j£j C Z"^, since in this case, the 
quantity considered is an integer. 

(iii) Consider the cubic case (7 = 1, and suppose that {Kj}j£j is included in a 
rectangular net. Up to translation, this rectangular net has the form 

{Am e I m e Z"^}, 
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with A a dxd matrix of the form A = RD, where D is diagonal, and Ris a rotation. 
Then we have, for all k,l,m £ U^: 

\Ak -AI + Amf - \Akf + \Alf - \ATaf = - Al) ■ {Ak - Am)\ 

Since ^ AA = D^, denoting /i^, . . . , /i^ the squares of the eigenvalues of D, Assump- 
tion 6.2 is then satisfied if and only if the group generated by fif, . . . in R is 
discrete, i.e. these numbers are (pairwise) rationally dependent. 

The reason for imposing the above assumption is a small divisor problem, as 
will become clear from the proof of the following lemma. It is possible to relax 
Assumption 6.2 to a less rigid one, to the cost of a more technical presentation. 
The latter is sketched in an appendix. 

Lemma 6.4. For = T'' or M ~ R'*, let rf be given by (6.3) and denote 

(6.5) Rl{t,x):=[ U%t~T)rl{T,x)dT, on[0,T]xM. 

Jo 

Let Assumptions 6.1-6.2 hold. Then, there exists a constant C > 0, such that: 

II-RiIIl°=([o,t];H'(7M)) ^ Ce. 

Proof. We only treat the case on = R"^ in detail. The case Ai ~ can be 
treated analogously. We have 

Rl{t,x) = 

E rC/^(t-T)((a,,a,,...a,,„^Je<^^i-^^=+-+'*^-+0/-) (r,x)dT. 

Thus, setting bi^^,,,_i^^^^ := ai-^ai^ . . .a^j^+i. Lemma 5.7 yields 

P?IIl~([0,T];V1/) £ £ TTi ~l 7(ll^<?i,...,Willi°°([0,T];Li) 

(fi,...,£2<, + i)6W 2'^+^^ 

+ l|A6£i,.../2„+J|ioo([o^T].Li) + ||9t6£i,...,^2,.+ilU°°([0,T];Li)) 
£e X! (ll^fi,---,«2^+illL=°([0,T];Li) 

(fi,...,£2<T+i)eJV 



where we have used Assumption 6.2. Next, using Young's inequality, as in the proof 
of Lemma 5.11, we get: 

X] ll^fi,...,«2„+il!L°=([0,T];Li) ^ X! ll"«ilU??Li • ■ • l|a«2^+illL~Li 

(«i,...,f2<,+i)6Af (^i,...,£2<, + i)eAr 



X! l|a<?ilUffLi • ■ • l|a<?2<,+ilU5?Li 



(fi,...,^2„+i)eJ2-+i 

^ ll(%)je-f|lL°^([0,T];_E)- 

Leibniz formula and Holder inequality yield similar estimates for A6£j^^...^£2<t+i ^-i^d 
3^^i,...,«2<,+i in i°°([0,r];ii(A^)), and the lemma follows. □ 
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6.2. Accuracy of the multiphase WKB approximation. With the above re- 
sults in hand, wc can now prove our main theorem. 

Theorem 6.5 (General approximation result). Let a ^ I, Ai = T"^ or M."^ , and 
Assumptions 6.1-6.2 hold. Given an approximate solution u|pp G C([0,T]; W{M)) 
as in (6.1), we consider a family of initial data (uq)£>o G W{A4), such that 



^0 '"app|t=0 



W{M) 



for some Cq ^ independent of e. Then there exists So{T) > 0, such that for 
any < £ ^ eo(T), the exact solution to the Cauchy problem (5.2) satisfies G 
i°°([0, T]; W^(A^)). In addition, u|pp approximates up to ©(e).' 

II, ,e „ ,,e II < |L,e _ ,,e || < np 

ir app||ioo([o,T]xA1) H app||ioo([o,T];M''(A1)) ' 

where C is independent of e. 

Obviously the result for a; € T, announced in the introduction, can be seen as a 
special case of Theorem 6.5. 

Proof. From Proposition 5.8, we may consider a solution G C([0, T^], W(A^)) 
to (1.6). We define the difference := - <pp. Then € C([0, t^], iy(7W)), 
where r"^ = min(T'^, T). We prove that for e sufficiently small, w'^ may be extended 
up to time T, with € C{[0,T],W{M)). Take Eq > so that CqEq < 1/2, and 
for e g]0, £o], let 

:=supft e [0,T] I sup \\w'it')\\w^M)^^}- 
t'e[o,t] 

We already know that > by the local existence result for it^. By possibly 
reducing Eq > 0, we shall show that T. The error term solves: 

idtw' + |Aw;^ = A (Kpp + «;1'"«pp + w') - |<ppP'<pp) + \r\ ~ erj, 

where rf , r| are given in (6.2)-(6.3). Using Duhamel's formula we can rewrite this 
equation as 



'{t) = W{t)wl ~ tX / W{t - t) (Kpp + i/;l2'^«pp + wn - |<ppr<pp) (r) dr 

JQ 

-iXRl{t)+ie [ W{t~T)rl[T)dT, 



where R\ is defined in (6.5). Using the fact that U^{t) is unitary on W{M.), and 
the estimates given in (6.4) and in Lemma 6.4, we obtain on [0, i^]: 

W{t)\\w(M) ^ Cie+|A| / |l(|<pp + u'1'"«pp + i«")-|<ppPXpp)(T)lk(A^)dT 



^Cie + Ca / ||w^(t)|1tv(7W) rfT, 

JO 

by the Lipschitz property from Lemma 5.5 . Note that, in view of Lemma 5.2 , resp. 
Lemma 5.4, (u|pp)e>o is a bounded family in C([0, T],W{M.)), and restricting t to 



24 



R. CARLES, E. DUMAS, AND C. SPARBER 



[0,t^] ensures that w'^{t) stays bounded in W{M). The constants Ci, C2 depend 
on Co and u|pp. Now, Gronwall lemma yields 



\\w'{t)\\w(M) ^ Cie (^1 + 



C2 



and we may reduce eg so that Ci^o (l + eP'^'^ jC-i) < 1. This shows that ^ T, 
for all e g]0, Eq]. Then, T"^ ^ T follows, as well as the desired approximation of u'^ 
by < since = 0{e) in L°"{[0,T];W). 



a 



7. Proof of the instability result 



This section is devoted to the proof of Theorem 1.2. To this end we essentially 
rewrite the proof of M. Christ, J. CoUiandcr, and T. Tao [n] in terms of weakly 
nonlinear geometric optics. It then becomes easy to sec that the justification given 
in the previous paragraph makes it possible to extend the one-dimensional analysis 
of [8] in order to infer Theorem 1.2. 

Proof of Theorem 1.2. We start with two Fourier modes, one of them being zero: 

idtu+^Au = X\u\^'^u ; u{0,x) ^ ao + aie''^''\ KeN. 

The fact that we privilege oscillations with respect to the first space variable is 
purely arbitrary. Define u as the solution to the same equation, with data 

u{0,x) = So + Sie'-^^^ 

Let 

1 ..fx 
e=— ; u {t,x)^ui^t,^ 

(e is chosen so that we remain on the torus.) We see that u'^ solves (1.6) on T"^, 
with 

From Theorem 6.5, we know that there exists T > independent of e, such that 

11"" - "IppIIl~([0,T]xT'') + 11^^ - ^|pp||i=c([o,T]xT'') = 

where M|pp is the approximate solution defined by (6.1), and u|pp is defined simi- 
larly. On the other hand, we have 



u {t, Kx) . 



a\e 



CT+1 

n 



aoe 



<pp(i,x) ^^aoe-'^*"" 
where, in view of (4.3), 9q is given by 

n=0 

We infer, uniformly in i G [0,T], 

/ {u{t,x) — u{t,x))dx 

J-fd 

with obvious notations. 

To prove the first point of Theorem 1.2. set 



~i\tei i(xi-t/2)/e 



|aoP""'"'|ai|'"- 



iXtBa 



aoe 



iXtBo 



0{e), 



~ P P 
ao = ao = -, «i = ^ 
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We infer, for < (5 < 1, 



> 



We have ||w(0) — u{0)\\h'^ < S provided K > 5^/*. Since uq = ao, we also have 



{u{t, x) — u{t, x))dx 



2ao sin 



0{e). 



We infer that we can find t e [0, S\ so that the right hand side is bounded from 
below by p/2, provided N is sufficiently large (hence £ sufficiently small). 

To prove the second point of Theorem 1.2, set 



ao 



ao = ao + d, ai ai = — — = - A ' 

IK * z 



For 5 small compared to p, we use the same estimate as above. 



{u{t, x) — u{t, x))dx 



> 



2ao sin 



for K sufficiently large. We now have 



70 — Oq 



> 



Pol 



l"ll 



>S+(pK 



2cr-2 



Now we see that if we assume a ^ 2, the left hand side can be estimated from below 
by 1/(5, provided N is sufficiently large, and we conclude like for the first point. 

To prove the last point in Theorem 1.2, we resume the argument of [2.^]. Fix 
ao e C \ {0}, and let ai G C to be fixed later. As K ^ oo, we have: 



u(0, •) ^ ao =:w(0,-) weakly in L2(T'i) ; ||it(0)||i2 ^ |ao|^ 
For any t > 0, we have, as K cx), 

u(t,x) aoe-'^*^° weakly in L^iT'^), 

where 



l"i| 



cr + 1 

n 



\ao\ \ai\ . 



Note that for any ao G C \ {0} and any angle 9 £ [0, 27r[, we can find ai G C so 



that 9o = e 
is given by 



We infer 



|ao 



2cr 



On the other hand, the solution to (1.8) with initial data ao 



u(t, x) = aoe' 



iXtli 



w- Van u{t,x) ~ u{t,x) = aoe-'^*^"°^ 



-ixte 



1 



For all t ^ 0, one can then choose 9 so that Xt9 ^ 27rZ. The discontinuity at ao of 
the map ao u{t), from L^{T'^) equipped with its weak topology into (C°°(T'^)) , 
follows. □ 
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Appendix A. A more general set of initial phases 
We can actually replace Assumption 6.2 with the following more general one: 
Assumption A.l. There exist 6 ^ 0, c > such that for all . . . , £2<t+i) G N, 



2a+l 



2a+l 
p=l 



satisfies: 



In §6, we have considered the case 6 = 0. However, allowing constants b > 0, 
we show that the assumption is satisfied by wave vector sets included in generic 
finitely generated nets. 



Proposition A. 2. For all p € W, there exist C,b > and Z C R ^ with zero 
Lebesgue measure such that, for all {ki, . . . , Kp) G M.'^^ \ Z , the set {Kj)j(zj con- 
structed from these initial wave vectors {Kj}j^jg satisfies Assumption A.l. 

Proof. We shall prove that the above result holds when Assumption A.l is replaced 
by the stronger one, where N is replaced by J'^"'^-^. 

All the wave vectors we consider belong to the group generated by {^i, . . . , Kp}. 
Thus, to each £k & J corresponds (afe.i, ■ • ■ , ctk.p) G such that: k^^. = ak,iKi + 
■ ■ ■ + ak,pKp. With this notation, for all {£i, . . . , £2^+1) G J^'^^^, we have: 



> £2<j+l) 



2a+l 



2a+l 



J2 + E (-ir+'lE"™^^'^^ 

j = l m=l j=l 



k=l 



2a+l 



m—1 
2a+l 



E E - E (-i)"«m,.«„ 

ij = l \k,e=l m=\ 

Now, a standard Diophantine result (see 6.17. [1, 11]) ensures that, for all choice 
of (Ki ■ Kj)i^ij<jp but in some subset of W with measure zero, we have, for some 
6' ^ and C''> 0: 

I p \ 

E ^i^i'^i ■ 

Such an estimate is then valid for almost all . . . , Kp) in (M'^)^. We apply it with 

2cr+l 2(T+1 



V(A,j) 



eZP \{0}, 




A. 



E(-i) 



k,e=i 



so that 



2cr+l 



E l/3.,.K2 E l«fe.-IK' 



k.i=l 



k=l 
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Now, choosing ki, . . . , Kp Q- linearly independent (which is true almost surely), we 
get that there exists a constant c > such that 

d 

Va e Q'', |ai| H h \ap\ < Kai^i H h Q!pKp)j|. 

Increasing c if necessary, so that c ^ 1, we get, when ki^ = ak^iKi + • • • + ak^pKp'. 
(a/c, ) ^ c{ki^). Finally, using the constants b' and C" from above, the desired 
estimate follows with b = 2b' and C = (2(2cr + ifc^y^'C. □ 

Under Assumption A.l (which is fairly general for plane waves, in view of the 
above proposition), we can easily adapt the analysis of §6. Essentially, we have to 
(possibly) strengthen the assumptions on the initial profile, in the case of = R"*, 
where we generalize Assumption 6.1 to: 

Assumption A. 3. On A4 = M.'^, the initial amplitudes satisfy: 
VI/3K2, ((K,)''9f«,)jeJo ei?(K'^), 

From Proposition 5.12, these data produce a solution {aj)j^j E C([0, T], E{Ai)) 
to the profile system (5.3). We consequently define the approximate solution it^pp 
as before 

<pp(^,a;) = ^a5(^,x)e*(*'-)/^ 
where the sequence {aj)j^j is now such that 

((^,)'afa,) ec([0,r],i?(x)), |/?K2, 

\ ^ je.J 

We can then reproduce the analysis of §6: Lemma 6.4 is still valid under Assump- 
tion A.l and A. 3, by straightforward verification. Then one just has to notice that 
this is the only step where the absence of small divisors plays a role in the proof 
of Theorem 6.5. Therefore, Theorem 6.5 remains valid under Assumption A.l and 
A.3. 
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